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Abstract 

Hot hypernuclear matter is investigated in an explicit SU{3) quark model 
based on a mean field description of nonoverlapping baryon bags bound by the 
self-consistent exchange of scalar a, C and vector uj, (p mesons. The a, u mean 
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fields are assumed to couple to the u, d-quarks while the £, <f> mean fields are 
coupled to the s-quark. The coupling constants of the mean fields with the 
quarks are assumed to satisfy SU(6) symmetry. The calculations take into 
account the medium dependence of the bag parameter on the scalar fields 
cr, We consider only the octet baryons N, A, E, H in hypernuclear matter. 
An ideal gas of the strange mesons K and K* is introduced to keep zero net 
strangeness density. Our results for symmetric hypernuclear matter show that 
^ \ a phase transition takes place at a critical temperature around 180 MeV in 

which the scalar mean fields cr, £ take nonzero values at zero baryon density. 
Furthermore, the bag contants of the baryons decrease significantly at and 
^ . above this critical temperature indicating the onset of quark deconfinement. 

The present results imply that the onset of quark deconfinement in SU{3) 
hypernuclear matter is much stronger than in SU(2) nuclear matter. 
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I. INTRODUCTION 



The majority of nuclear phenomena have been successfully described in relativistic mean- 
field theory using only hadronic degrees of freedom • However, due to the observations 
which revealed the medium modification of the internal structure of the baryons ||, ft 
has become essential to explicitly incorporate the quark-gluon degrees of freedom while 
respecting the established model based on the hadronic degrees of freedom in nuclei. One 
of the first models put forward along these lines is the quark-meson coupling (QMC) model, 
proposed by Guichon ||, which describes nuclear matter as a collection of non-overlapping 
MIT bags interacting through the self-consistent exchange of scalar a and vector oj mesons 
in the mean field approximation with the meson fields directly coupled to the quarks ||[7|. 
The scalar a meson is supposed to simulate the exchange of correlated pairs of pions and 
may represent a very broad resonance observed in tttc scattering, while the vector uj meson 
is identified with the actual meson having a mass of about 780 MeV. In the chiral models, 
the scalar a and vetcor to mean fields represent the u, d quark condensates @|J. The QMC 
model thus incorporates explicitly the quark degrees of freedom and this has nontrivial 
consequences ||[7]|. It also have been extended to study superheavy finite nuclei [jnjfnj 



In the so-called modified quark meson coupling model (MQMC) it nas been 



further suggested that including a medium- dependent bag parameter may be essential for 
the success of relativistic nuclear phenomenology. It was found that when the bag parameter 
is significantly reduced in the nuclear medium with respect to its free-space value, large 
cancelling isoscalar Lorentz scalar and vector potentials for the nucleon in nuclear matter 
emerge naturally. Such potentials are comparable to those suggested by relativistic nuclear 
phenomenology and finite density QCD sum rules [14]. The density-dependence of the bag 
parameter is introduced by coupling it to the scalar meson field [[Tj|. This coupling is 



motivated by invoking the nontopological soliton model for the nucleon In this model 
a scalar soliton field provides the confinement of the quarks. This effect of the soliton field 
is, roughly speaking, mimicked by the introduction of the bag parameter in the Bag Model. 
When a nucleon soliton is placed in a nuclear environment, the scalar soliton field interacts 
with the scalar mean field. It is thus reasonable to couple the bag parameter to the scalar 
mean fields OJil 



The QMC model was extended to finite temperatures to investigate the liquid-gas phase 
transition in nuclear matter [16|. Recently, the MQMC model has been also extended to 
finite temperature an d nas been applied to the study of the properties of nuclear 

matter where it was found that the bag parameter decreases appreciably above a critical 
temperature T c « 200 MeV indicating the onset of quark deconfinement fT8"l . The effect 



of glueball exchange as well as a realization of the broken scale invariance of quantum 
chromodynamics has also been investigated in the MQMC model through the introduction 
of a dilaton field |T^[. It was found that the introduction of the dilaton potential improves 



the shape of the saturation curve at T=0 and affects hot nuclear matter significantly |1S 



In the present work, we extend the MQMC model to hot hypernuclear matter by intro- 
ducing the scalar ( and vector mean fields that are coupled to the s-quark in addition to 
the a and uj fields which are coupled to the u, d-quarks. The ( and <fr fields are identified with 
the real mesons having masses = 975 and = 1020 MeV respectively P,|9j|2"0]j . Hy- 
pernuclear matter is considered to contain the octet p, n, A, E + , S°, S~, H° and H~ baryons 
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pTfl . We introduce an ideal gas of kaons to keep zero net strangeness density ps = 0. We 



simplify the calculations by considering symmetric hypernuclear matter whereby the octet 
baryons reduce to 4 species: 2N, A, 3E, 2H. 

The outline of the paper is as follows. In section II, we present the MQMC model for 
hypernuclear matter at finite temperature, together with the details of the self-consistency 
conditions for the vector and scalar mean fields. In section III, we discuss our results and 
present our conclusions. 



II. THE QUARK MESON COUPLING MODEL FOR HYPERNUCLEAR 

MATTER 

The quark field ipqif^t) inside a bag of radius Ri representing a baryon of species i 
satisfies the Dirac equation 

^(r,t) = 0. (1) 

where m q is the current mass of the quarks of flavor q and where r refers to the up or down 
quarks and s refers to the strange quark. The Kronecker deltas insure that the u, d quarks 
are coupled only to the a and u fields while the s-quark is coupled only to the £ and fields. 
In the mean field approximation all the meson fields are treated classically and, for nuclear 
matter in this approximation, these fields are translationally invariant. Moreover, because 
of rotational invariance, the space-like components of the vector fields vanish so that we 
have ujfiY 1 * =< ojq > 7 = wy° and 4>^y^ =< <fio > 7 = 07°. 

The single-particle quark and antiquark energies in units of for quark flavor q are 
given by 

e nn = Q n K ± (g^R^ + g^R^) (2) 

where 

Q n q K = sJx q nK 2 + R\m*\ (3) 

and 

m* q = m° q - g q a a5 qr - g^ q „ (4) 

are the effective quark kinetic energy and effective quark mass, respectively. The boundary 
condition for each quark of flavor q at the bag surface is given by 

z 7 -nVC(*L)=< K «J, (5) 

which determines the quark momentum x q nK in the state characterized by specific values of 
n and k. For a given value of the bag radius Ri for baryon species % and the scalar fields a 
and (, the quark momentum x q nn is determined by the boundary condition Eq. (^j) which, 
for quarks of flavor q in a spherical bag, reduces to jo(^ K ) = Ajii(^ K ), where 

fi*™07,C) + ^m*(<T,C)' 1 ' 



i^d, -m° q + (g'a - gl^)8 qr + (g^ ~ gfa-ftf, 
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The quark chemical potential fi q , assuming that there are three quarks in the baryon 
bag, is determined from 3 = J2 q n q where n q is the number of quarks of flavor q and is 
determined by [p|-|T9| 



n„ 



E 



i 



i 



,(e q ™/R-» q )/T 



+ 1 



+ 1 



(7) 



The total energy from the quarks and antiquarks of each baryon of species i is 

1 1 



E 



tot 



E 

nnq 



n, 



Qnft r 

' R 



+ 



(8) 



The bag energy for baryon species i is given by 

El„„ = E, 



(9) 



where Bi = B^a, () is the bag parameter. In the simple QMC model, the bag parameter B 
is taken as B corresponding to its value for a free baryon. The medium effects are taken 
into account in the MQMC model [0,0 by coupling the bag parameter to the scalar meson 
fields. In the present work we generalize the coupling suggested in the latter references to 
the case of two scalar meson fields by using the following ansatz for the bag parameter 



Bi = B Q exp 



4:((n u + n d )g^a + n s gf() 
3 Mi 



(10) 



with gf? and as additional parameters. The spurious center-of-mass momentum in the 
bag is subtracted to obtain the effective baryon mass p2 
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and 



<x 2 > l =E' 
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e {e q ™/Ri-n q )/T + 1 1 e (e q ™/Ri+n q )/T + j 



(11) 



(12) 



(13) 



The bag radius Ri for baryon species i is obtained through the minimization of the baryon 
mass with respect to the bag radius |7|,p^,p!7HT9 



dM* _ Q 



dRi 



(14) 



The total energy density at finite temperature T and at finite baryon density pb reads 
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Baryons 1111 



(27T) 

+ eg, (15) 

where 7$ is the spin-isospin degeneracy factor of baryon species i and where the last term 
corresponds to the energy density of the if-mesons treated here as an ideal gas. In Eq. 
(|15|) fi and f i are the Fermi-Dirac distribution functions for the baryons and antibaryons of 
species i, 

fi = e «-^)/T + 1 ' ( 16 ) 

and 

7i = e «+M*)/T + p ( 17 ) 

where e* and //* are, respectively, the effective energy and effective chemical potential of 



baryon species i. These are given by e* = J k 2 + M* 2 and //* = B^b + Si/As 1 — (^o; + g$i4>) 
where Bi and Si are the baryon and strangeness quantum numbers and where = 
{n u + ridjglji an d 9<t>% = n s9\% ■ The chemical potentials [is, Us are determined by the 
self-consistency equations for the total baryonic density 

PB = J^yE B ^i J -7<), (18) 

and the total strangeness density 

Baryons Kaons 

(27T)" 



1 naryons Maoris 

/'>- = 7S3s E # 7 i/ -7 4 ) + E $p& = (19) 



where in the last equation we have introduced the contribution of an ideal gas of if and if* 
mesons to make the total strangeness density vanish identically. The vector mean fields are 
determined by 

gu 



<* = E ^B iPi , (20) 

and 

g<t> 



E^^A- (21) 



n? 



The pressure is the negative of the grand thermodynamic potential density and is given 

by 

P=-Y -1L- [ (P&fi + JA + -mlco 2 + -m> 2 - -mla 2 - -m 2 c ( 2 + P%, (22) 
3 Y (2?r) 3 i e* JtJ 2 2 2 2 c v ; 

where the summation i runs over the 8 species of the baryon octet which reduces for symmet- 
ric hypernuclear matter to 4 species with 7^ = 4, 2, 6 and 4 for N, A, S and 2, respectively. 
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In Eq. (p2[) is the pressure of the ideal gas of K mesons. The density of the .fT-mesons 
of species % is given by 

PKi = ^\ d'Kh - h) (23) 

where the spin-isospin degeneracy = 4,6 for K,K*, respectively. The Bose-Einstein 
distribution functions for the K mesons are given by 

h = — ; (24) 

and 

(25) 
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where fii = S^s is the chemical potential of i-T-meson of species i. The total energy density 
and total pressure of the K-meson ideal gas are given by 



e K = E J^js J <*V fc2 + M ^h - h), (26) 



K 



and 



P K=\T,l^yJ d 3 k f (h - 1,), (27) 
3 i (2vr) d J Jk 2 + M 2 



respectively. 

The scalar mean fields a and ( are determined through the minimization of the ther- 
modynamic potential or the maximizing of the pressure with respect to these fields. The 
pressure depends explicitly on the scalar mean fields through the \m 2 a a 2 and ^m 2 ( 2 terms 
in Eq. (p2|). It also depends on the baryon effective masses M* which in turn also depend 
on cr and (. If we write the pressure as a function of M* , a and ( [[H|,[n|, the extremization 
of P(M*, a, () with respect to the scalar mean field a can be written as 




where 




= -ml<r, (29) 

{M*} 

with a similar experssion for the extremization of P(M*,a,() with respect to the scalar 
mean field (. The derivative of the pressure with respect to effective mass M* reads 
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(30) 



Since the baryon chemical potential \xb and temperature are treated as input parameters, 
the variation of the vector mean field uj with respect to the effective baryon mass M* at a 
given value of the baryon density pB reads 



duo 
dMf 



[Qui/ 


m 2 ] [ ll /(2n) 3 ] [l/T]Jd 3 k^ 

i 


f i (l-f i )-f.(l-f.j 




l + E, 


glj/ m l_ 


[ Ti /(27T)3] [l/T\Jffik 


fiO--fi) + fiO--fi)_ 



(31) 



with a similar expression for {d<p/ dM*) T . The coupling of the scalar mean fields a and 
( with the quarks in the non-overlapping MIT bags through the solution of the point like 
Dirac equation should satisfy the self-consistency condition. These constraints are essential 
to obtain the correct solution of the scalar mean fields a and (. 



III. RESULTS AND DISCUSSIONS 

We have studied hypernuclear matter at finite temperature using the modified quark 
meson coupling model which takes the medium-dependence of the bag into account. We 
choose a direct coupling of the bag parameter to the scalar mean fields a and ( in the form 
given in Eq. fllPP - The bag parameter is taken as that adopted by Jin and Jennings [12| 
B^ 4 = 188.1 MeV and the free nucleon bag radius Rq = 0.60 fm. We have taken the current 
quark masses to be m u = m d = and m s = 150 MeV. For g q a = 1, the values of the vector 
meson coupling constant and the parameter g^, as fitted from the saturation properties of 
nuclear matter, are given as # 2 /47r = (3^) 2 /47r = 5.24 and gf' ''/Att = (3^ 9 ) 2 /47r=3.69. 
The Zi = 2.03, 1.814, 1.629 and 1.505 are chosen to reproduce the baryon masses at their 
experimental values Mat,a,s,h =939, 1157, 1193, and 1313 MeV respectively. Normal nuclear 
matter saturation density is taken as pbq = 0.17 fm -3 . The extra coupling constants needed 
to couple the scalar and vector mean fields ( and <fi to the s-quark are chosen to satisfy 577(6) 
symmetry where \g*\ = \f2\gl\ and \gV\ = y/2\g^\ and \g^ q \ = V2\g^ q \. If it is assumed that 
the mean fields ( and <fi are positive definite, then all the coupling constants are positive 
and the absolute value signs become redundant. 
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The a mean field is supposed to simulate the exchange of correlated pairs of pions and 
may represent a very broad resonances observed in nir scattering and fixed at m a = 550 MeV, 
while the vector uo meson is identified with the actual meson whose mass is m u = 783 MeV. 
Since the mean fields, a and u, are considered as < ud > condensates, they interact only 
with u, d-quark in the baryons. On the other hand, the scalar and vector mean fields (, <fi are 
considered as actual mesons with = 975 MeV and = 1020 MeV, respectively. They 
are considered as < ss > condensates and interact only with the s-quarks in the baryons. 
This picture is consistent with the chiral models [|,[| . We fixed the total strangeness density 
Ps to by introducing an ideal gas of K and K* mesons where rrix = 495 MeV and ttik* = 
892 MeV. The contribution of other K mesons was found to be negligible. It is supposed 
in the ideal gas limit that the Kaons do not interact with the mean fields ( and <fi. The 
extension to the case that the Kaons interact with the ( and <fi fields will be considered in 
a future work. 

We first solve Eqs. (p"8"D, (p"9"D, (|20|) and (|21|) self consistently for given values of tem- 



perature T and densities ps and ps to determine the baryonic and strangeness chemical 
potentials \ib and /is, respectively. These constraints are given in terms of the effective 
baryon masses M* which depend on the bag radii Ri, the quark chemical potentials p, l q 
and the mean fields. For given values of the scalar fields a, ( and vector fields lu, <f), the 
quark chemical potential and bag radius of species i are obtained using the self consistency 
conditions Eqs. (0) and ([[4]), respectively. The pressure is evaluated for specific values of 
temperature T and chemical potentials (/ie,/is) which now become input parameters. We 
then determine the values of a and ( by using the extermization conditions as given in Eq. 
(f28|). These constraints take into account the coupling of the quark with the scalar mean 
fields in the frame of the point like Dirac equation exactly ||7,T2| . 

The dependence of the baryon effective masses ash 011 the total baryonic density 
Pb and temperature is shown in Fig. 1 where it is seen that the baryon masses decrease 
with baryonic density except at the highest temperatures where the effective masses become 
almost density- independent. Moreover, for a given baryonic density ps it is seen that as the 
temperature is increased the mass M* of species i first increases slightly up to about T = 150 
MeV and then decreases rather rapidly for T > 150 MeV. These results are displayed in a 
different manner in Fig. 2 which plots the baryon masses M* as a function of T for ps = 
fm ~ 3 and pbq = 0.17 fm ~ 3 . It is seen that the effective baryon masses Mjy- A)SH with 
Pbq = 0.17 fm ~ 3 are less than those with zero baryonic density. Moreover, the effective 
baryonic masses M^ ATiE increase only slightly, if at all, with temperature up to about 
T = 150 MeV beyond which they decrease rapidly. This behaviour is qualitatively similar 
to our earlier results for normal nuclear matter fl^JTMHI where the rapid decrease in the 



nucleon's effective mass was, however, found to start at rather higher temperatures T > 200 
MeV. This rapid decrease of M* with increasing temperature resembles a phase transition 
at high temperatures and low density, when the system becomes a dilute gas of baryons in 
a sea of baryon-antibaryon pairs |23| . 

In Fig. 3, we display the scalar mean fields a and ( as functions of the total baryonic 
density ps for various temperatures. It is seen that the value of a initially decreases with 
increasing temperature for temperatures less than 150 MeV. The scalar mean field ( is 
almost negligible for such low temperatures. However, as the temperature reaches 150 MeV 
there are indications of an increase in a at low baryon densities where it attains a nonzero 
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value at ps = 0. For still higher temperatures, the situation is more dramatic with the 
value of a increasing with temperature for all values of pB- This is also qualitatively similar 
to our earlier results for SU{2) nuclear matter [|1^,|T^]. Furthermore, the scalar mean field 
( becomes important for T > 150 MeV and also increases with temperature for all values 
of pb- An interesting new feature here is that the scalar mean fields a and ( tend to take 
almost constant values irrespective of density at high temperatures which was not seen in 
our earlier calculations for SU(2) nuclear matter even at temperatures as high 240 MeV . 
Fig. 4 displays o and ( versus T for p B = and pb = 0.17 fm ~ 3 . It is seen that the a field 
has a nonzero (and almost constant) value only at the higher density until the temperature 
reaches about T = 150 MeV when a rapid increase sets in at both densities. The increase at 
zero density is actually more dramatic and the a field rapidly attains values equal to those 
occuring at ps = 0.17 fm ~ 3 . This is another indication of a phase transition to a system 
of baryon-antibaryon pairs. The behaviour of the ( field is qualitatively similar except that 
its value is negligible at both densities for temperatures less than about T = 150 MeV. 

In Fig. 5, we display the baryonic density dependence of the bag parameters for N, A, S, H 
for different values of the temperature. For each baryon, the bag parameter increases with 
temperature for temperatures less than 150 MeV. However, the situation is completely re- 
versed after the phase transition takes place. For temperatures T > 150 MeV the bag 
parameters display a dramatic decrease with temperature for all densities. This can be seen 
more clearly in Fig. 6 which displays Bi vs T for pb = fm ~ 3 and p# = 0.17 fm ~ 3 . 
The bag parameters are almost constant until the temperature exceeds T = 150 MeV when 
they start to decrease rapidly. This indicates the onset of quark deconfinement above the 
critical temperature: at high enough temperature and/or baryon density there is a phase 
transition from the baryon-meson phase to the quark-gluon phase. This behaviour is also in 
qualitative agreement with our earlier results for 577(2) nuclear matter I?| except that the 
decrease here is more dramatic indicating that the phase transition is much stronger than 
in ordinary nuclear matter. Our results are also comparable to those obtained from lattice 
QCD calculations which have so far only explored the zero baryon density axis of the phase 
diagram in a meaningful way. The lattice results with 2 light quark flavors, indicate that the 
transition from hadronic matter to a quark-gluon plasma occurs at a temperature T « 140 
MeV and that for low densities it may not be a phase transition at all but what is called a 
rapid crossover p^j . 

Fig. 7 displays the relative abundance Pi/pB for each baryon species. At low temper- 
atures the nucleons N are almost the only constituents. However, the contribution of the 
hyperons starts to be noticeable when the temperature reaches 100 MeV and becomes more 
important when the temperature is increased to T > 150 MeV. At temperatures T > 200 
MeV the contribution of the nucleons falls down to about half the total baryonic density. 
This can also be seen in Fig. 8 where we display the ratio of the baryonic strangeness den- 
sity to the total baryonic density. As mentioned earlier, we keep the net strangeness of the 
system fixed to zero by introducing Kaons so that ps = p^ aryons -\- p^ esons = o. It is seen that 
the net strangeness of the baryons p^ aryons [ s small at low temperatures T < 150 MeV (note 
the logarithmic scale). However, as the temperature increases the net strangeness of the 
baryon octet becomes significant. For T > 150 MeV, the ratio p^ aryons / Pb becomes of order 
one. Therefore, the .fT-meson plays a significant role at high temperatures and probably also 
in the phase transition. 
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In Fig. 9, we display the relative abundances of the various baryon species p%/ Pb versus 
temperature at normal nuclear matter density p# = 0.17 fm ~ 3 . It is seen that the the 
nucleons are dominant at low temperatures and their relative abundance Pn/pb decreases 
very slowly at first and deviates very little from 1 until the temperature reaches about 100 
MeV when pn/pb starts to decrease rapidly and becomes negligibly small for temperatures 
larger than about 280 MeV. On the other hand, the abundance of the hyperons is negligible at 
low temperatures and increases significantly as the temperature is increased beyond T = 100 
MeV. They become more abundant than the nucleons for temperatures larger than 200 MeV. 
The relative abundance of hyperons in high energy heavy ion collisions can therefore be used 
as a simple thermometer to measure the temperature of the hot nuclear matter produced in 
the reaction and probably to study the phase transition to the quark gluon plasma. 

In conclusion, we have generalized the MQMC model to the case of hot hypernuclear 
matter by introducing two new meson fields that couple to the strange quarks and using 
SU (6) symmetry to relate the coupling constants. The results are qualitatively similar to 
those obtained for SU(2) nuclear matter including the onset of quark deconfinement. It 
is observed, however, that in this model quark deconfinement in SU(3) (or hypernuclear) 
matter is much stronger than in 577(2) (or normal) nuclear matter and that hyperons become 
more abundant than nucleons at high temperatures T > 200 MeV. An ideal gas of Kaons 
was introduced to keep zero net strangeness density. The interaction of the Kaons with the 
meson fields will be considered in a future work. 
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FIGURES 



Fig. 1 





FIG. 1. The effective baryon masses M* in hypernuclear matter as a function of the total 
baryonic density pg for different values of temperature. 



12 



Fig. 2 




FIG. 2. The effective baryon masses M* in hypernuclear matter as a function of temperature 
for p B = and p B = 0.17 fm ~ 3 . 
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Fig. 3 




FIG. 3. The mean scalar fields a, ( as a function of the total baryonic density ps for different 
values of temperatures. 
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Fig. 4 




FIG. 4. The mean scalar fields a, ( versus the temperature for ps = and ps = 0.17 fm 
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Fig- 5 





FIG. 5. The bag parameters Bi of the nucleoli and hyperons versus the total baryonic density 
for different values of temperature. 
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Fig. 6 
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FIG. 6. The bag parameters Bi of the nucleoli and hyperons versus the temperature T for 
p B = and p B = 0.17 fm ~ 3 . 
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Fig. 7 
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FIG. 7. The relative abundance of baryon species i, j^, versus the total baryonic density 
for different temperature. 



18 



Fig. 8 




FIG. 8. The absolute value of the ratio of the baryonic stangeness density to the total baryonic 
density ps/pB as a function of the total baryonic density pb for various temperatures. 
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Fig. 9 
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FIG. 9. The relative abundance of baryon species i, j^, versus temperature T for ps = 0.17 
fm ~ 3 . 
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